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ABSTRACT 



Black hole accretion disks have been proposed as good candidates for a range of interesting 
nucleosynthesis, including the r-process. The presence of the black hole influences the neutrino 
fluxes and affects the nucleosynthesis resulting from the interaction of the emitted neutrinos and 
hot outflowing material ejected from the disk. We study the impact of general relativistic effects 
on the neutrinos emitted from black hole accretion disks. We present abundances obtained by 
considering null geodesies and energy shifts for two different disk models. We find that both the 
bending of the neutrino trajectories and the energy shifts have important consequences for the 
nucleosynthetic outcome. 

Subject headings: accretion disks, nucleosynthesis, neutrino physics, null geodesies 



1. Introduction 

Rapidly accreting disks around black holes arise 
in several scenarios. Disks with MeV or higher 
temperatures have been discussed in the context of 
compact object mergers, long and short duration 
gamma ray bursts, and core collapse supernovae 
which have a rotating progenitor. 

The coalescence of compact objects produces 



susceptible to detection with interferometric grav- 



itational wave detectors, e.g. Abbott (2008) and 
references therein. Energy release from these ob- 
jects has been suggested as a promising source for 



the generation of gamma ray bursts (Popham et 



aL][T999l [Ruffert fc Jailka|p99| [Setiawan et al 



hot disks ( 


Kluzniak & Lee 


1999 


Rosswog et al. 


2004 


Taniguchi et al. 


200£ 


), and the evolution 



of these objects is of great importance for sev- 
eral reasons. Mergers generate gravitational waves 



2006[|Kneller fc McLaughlin |2006[ ). The neutrino 
flux from these objects is so large, that it would 
be easily detected by currently online neutrino de- 
tectors ( McLaughlin&^^Surman[20^ Caballero et 



lcaballe@uw.edu. olca- 
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al]|2009 ) 

These objects will also eject nucleosynthetic 
products that must be considered when study- 
ing the galactic inventory of elements. Thus it 
is important to determine the type of elements 
that are produced from accretion disks. Collisions 
have been speculated to contribute to the synthe- 
sis of neutron rich nuclei via the r-process, due to 
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the decompression of cold or mildly heated neu- 



tron star matter during the merger (Lattimer & 



Schramm] 1 1976} |Meyer||1989[ [Goriely et al.||2005 



Freiburghaus et al.||1999[ [Metzger et al.||2010[ ). Tt 

has also been suggeted that they form r-process 
nuclei from the ejection of hot material during the 



merger (McLaughlin & Surman 2005; Oechslin et 
ar]|2007[ [Surman et al 112008; |Metzger et al.pOOg 



Wanajo fc Janka||2010[|Kajino et al.||2010[ ) 



McLaughlin & Surman (2005) and Surman et 



al. ( 2008 ) found neutron rich material in the out- 
flow ejected out of the plane of the inner hot re- 
gions of an accretion disk. A successful r-process 
occurred when either (1) low entropy and fast out- 
flow conditions obtained or (2) a favorable combi- 
nation of electron neutrino and anti-neutrino spec- 
tra, which when taken together with the accre- 
tion disk geometry allowed the material to flow 
far away from the source of the neutrinos before 
nuclei begin to form. Disks can form other ele- 
ments as well. Disks that are expected to occur 
in rare core collapse supernovae. e.g. "collapsars" 



(MacFadyen & Woosley 1999) are dominated by 



trapped electron neutrinos or are mostly trans- 
parent. This leads to elements that are formed 
when roughly equal numbers of neutrons and pro- 
tons are present, such as Nickel-56 and p-process 



elements ( 


Surman et al. 


2006 Pruet et al. 


2003 


Kizivat et al.||2010 


). 



The presence of a massive object such as a 
black hole changes the properties of space-time 
around it affecting the spectra of radiation emerg- 
ing from the matter located in its vicinity. Since 
the neutrino spectra are crucial to the outcome 



of the nucleosynthesis in hot outflows ( 


Meyer 


1995 McLaughlin & Fuller 1996| McLaug 


ilin et 


al. 1996 Haxton et al. 1997 Meyer et al. 1998 


Frohlich et al. 


2006 


) , particularly in accretion disk 


hot outflows ( 


Surman & McLaughlin 


2004 


), it is 



important to consider the effects of general rela- 
tivity on their energies and trajectories. 

Several detailed studies have been conducted 
for photons that emerge from accretion disks and 
these results are also applicable to neutrinos. In 



an early work Luminet (1979) studied the pho- 



ton spectrum from black hole accretion disks as 
seen by an observer located at infinity. The ef- 
fects of light bending and energy shifts for photons 
emerging from a Schwarzschild black hole were in- 



and observers located at different inclination an- 
gles with respect to the axis of symmetry of the 



disk. Fukue & Yokoyama ( 1998 1 studied the spec- 



tra of a non-monochromatic distribution by plac- 
ing the observer at an arbitrary distance from the 



black hole. Cunningham ( 1975 ) calculated the X- 



ray spectrum emerging from a disk around a Kerr 
black hole. Bhattacharyya et al. (20011 studied 



the relativistic spectra when the black hole is re- 
placed by a rapidly rotating neutron star. A more 
recent work compared different models for the ac- 
cretion disk and applied their results to ultralumi- 



nous X-ray sources (Lorenzin & Zampieri 20091. 
To date, studies of the effects of general relativ- 
ity, including ray bending, on neutrinos emerging 
from accretion disks have focused on neutrino pair 
annihilation rates and their production of gamma 



ray bursts (Asano & Fukuyama 2001 
2^|Zalamea fc Belobodorov||2010| ). 



Birk et al. 



eluded. Luminet (1979) also considered flat disks 



We consider neutrino general relativistic effects 
with a focus on the consequences for the produc- 
tion of nuclei in the vicinity of accretion disks. 
We are interested on the nucleosynthesis occur- 
ring when hot winds of free nucleons ejected from 
the disk interact with the neutrinos emitted by 
the disk. In this paper we incorporate effects such 
as neutrino bending and energy shifts to the neu- 
trino spectra. Our study includes two accretion 
disk models: a dynamical model coming from 3D 
simulations and a steady state disk in one dimen- 
sion. 

This paper is organized as follows: in Sec. [2] 
we present the basic formalism used to determine 
neutrino ray bending, energy shifts and neutrino 
fluxes. In Sec. [3]we introduce the disk models and 
describe our outflow model. In Sec. |4]we discuss 
the nucleosynthetic outcomes for various scenarios 
and in Sec. [5] we conclude. 

2. General Relativistic effects 

Our treatment of general relativistic effects is 
divided in two main components: bending of neu- 
trino trajectories and energy shifts. For the last 
effect we have found energy shifts in the Kerr and 
Schwarzschild metrics. However, for simplicity, we 
use only the Schwarzschild metric for our treat- 
ment of the neutrino trajectories. In the next two 
subsections we describe our approaches to each of 
these effects. We label the points from where neu- 
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trinos are emitted as rem and the points where the 
neutrinos fluxes are observed as rot- The emission 
points, Tern, correspond points on the neutrino sur- 



faces (found in a previous work (Caballero et al. 
2009 )) and rob to all the points on the outflow tra- 



jectory where nuclear products are synthesized. 
2.1. Neutrino trajectories 



Neutrino trajectories correspond to null geodesies 
in the curved space. Here we use the same nota- 
tion and follow a similar methodology as in |Muller| 
(2009) to find neutrino trajectories. Our main ob- 



jective is to find the angles ^ and a that a neutrino 
forms on the sky of an observer located at a dis- 
tance from a black hole. The angle ^ describes 
the direction of momentum of the neutrino with 
respect to line joining the observer and the black 
hole, while a is the angle with respect to a hor- 
izontal, both directions defined in a coordinate 
system centered at the observer. The angles ^ and 
a are found after solving the equation of motion 
for a neutrino in the vicinity of a black hole. 

We start by describing the neutrino trajectory 
in a coordinate frame centered at the black hole. 
Then we find the required variables in the ob- 
server's coordinate system. For simplicity we as- 
sume that the space time around the black hole is 
described by the Schwarzschild metric, 



-{1- — ]dt^ 
r 



{dO 



2 ' sin^edip^) 



1 

1 - rs/r 

2 



dr^ 



(1) 



where = 2M and M is the black hole mass. The 
spherical symmetry of the Schwarzschild metric 
makes it possible to solve the equations of motion 
in the = ti/2 plane. In this way the Lagrangian 
formalism results in null geodesies described by 



l_ fdr_\ 
\dLp J 



1 , , X 1 

+ 12 



62^ 



or 



dip = ± 



dr 



(2) 



(3) 



The sign of equation [3] depends on whether the 
neutrino is approaching or leaving the black hole. 
In the above equation 6, the impact parameter at 
infinity, is the ratio between the neutrino energy E 



and the its angular momentum L, b — L/E. The 
impact parameter 5 is a constant along the neu- 
trino trajectory and makes it possible to calculate 
the angle ^ between the neutrino trajectory and 
the radial direction at any point r by 



r sin^ 
a/1 - rs/r 



(4) 



A neutrino moves radially if 6 = and tangentially 
if ^ = 7r/2. At the point of closest approach to the 
black hole 



b = 



(5) 



At r+, the largest root of the square root in 
equation [3] there is inversion of the radial mo- 
tion. It corresponds to instantaneous tangential 
motion. Also, there exits a critical impact param- 
eter bcrit — V^Tcrit, with Tcrit = 3rs/2 such that 
null geodesies with b < bcrit are captured by the 
black hole. 

We can use equation |3] to our convenience. For 
example, we can start from an initial emission 
point (rem, Vem) and a known impact parameter b 
and draw neutrino trajectories starting from that 
point. Figure [T] shows null geodesies followed by 
neutrinos leaving at ipem — 0° and different dis- 
tances rem = 1^+ from a black hole. The impact 
parameter b is calculated using equation [5j 

We can also fix the initial and final points of the 
trajectory, which are the limits of the integrals in 
equation [3j and find the impact parameter b that 
satisfies it. By solving for b we are also finding the 
angle ^ that the null geodesic makes at any point 
of the trajectory by using equation [4] Depending 
on the position of the emitter and the observer 
with respect to the black hole this could imply 
solving consistently equation [5] to obtain the point 
of closest approach. Figures [2] and [3] show different 
trajectories for which we have found ^ using the 
procedure described above. 

In the present study the distance from the emis- 
sion points on the neutrino surface to the black 
hole are such that r^ = r^m > "fcrit- Therefore 
we are not concerned about geodesies that end up 
into the black hole. We are also interested in null 
geodesies that start from a fixed point on the neu- 
trino surface and end at another fixed point on 
the outflow. In such case we focus on flnding b 



3 



-10 



-5 
X [km] 



5 



X [km] 



Fig. 1. — Null geodesies, found by using equation 
[3j starting at Lpem = 0° and distances x = rem = 
r+ = 2.5,2.7,3,4,5,6 km. The black hole is lo- 
cated at {x = 0,y = 0) and its mass is 0.5Mq. 
The impact parameter b is found using equation 

El 

from equation |3] depending on the initial condi- 
tions rerm^em and rob,fob- Taking into account 
the relationship between the distances rem and rob, 
our study requires the solution of the following 
cases, which are sketched in Figure |4] (a study of 
more cases in the strong deflection limit can be 
found in Bozza & Scarpetta ( 2009[ )): 

Case 1. The distance from the observer to the 
black hole is larger than the distance from the 
emission point to the black hole, rob > ^em, and 
the neutrino does not reach the point of closest 
approach, r^m >'''+■ In this case neutrinos are 
leaving the black hole and we need to solve equa- 
tion [3] with the positive sign. 

Case 2. The distance from the emission point 
to the black hole is larger than the distance from 
the observation point to the black hole, r^m > fob, 
and the neutrino does not reach the point of closest 
approach to the black hole, rob > In this case 
the neutrinos are approaching the black hole, so we 
need to solve equation |3] with the negative sign. 



dip ~ 



dr 



b^ 



(6) 



Case 3. The distance from the emission point to 
the black hole is larger than the distance from the 
observer to the black hole and the neutrino trajec- 



Fig. 2. — Null geodesies, found by using equa- 
tion [3| starting from the vertical line at different 
distances r from a black hole (centered dot) and 
ending at the same observation point rob — 10 km, 
iPob = 0. The black hole mass is O.5-M0. ^ is the 
angle formed line joining the black hole and the 
observer and the direction of the null geodesic at 

Tob- 



tory reaches the point of closest approach to the 
black hole. In this case we need to split equation 
[3] in two branches. One is the approaching phase 
(neutrinos leave the neutrino surface and approach 
to r_|_), and the departure phase in which neutri- 
nos leave r_|. to get to the observer. We assign the 
signs in equation [3] accordingly to get 



dip = — 



dr 



62 



dr 



(7) 



This last case is similar to the situation where 
fern < rob and r+ is reached. For this latter situ- 
ation, one would switch rob and r^m for case 3 in 
Figure |4] 

Whether the neutrino reaches r^ or not is not 
known a priori. Therefore we need to consider 
all the cases when finding null geodesies for our 
problem. We use a root finding algorithm to solve 
the corresponding cases of equation |3] We find 
values of b and r+ such that equation [5] is valid 
and consistent with the appropriate case. 
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Fig. 3. — Null geodesies starting at different angles 
but from the same distances — 3rs to the 
black hole (BH) located at (x = ?/ = 0). The 
black hole mass is O.5M0 and = 2.9 km. The 
observer is at rob^lOr^. 




BH 



Fig. 4. — Different relationships between the dis- 
tances from the emission point rem to the black 
hole (BH) and from the observation point rot to 
the BH. r-|_ represents the point of closest ap- 
proach to the BH.The labels (1), (2), and (3) rep- 
resent the cases discussed in the text. 

Now we proceed to determine the angles that 
a neutrino describes on the observer's sky. Our 
purpose is to find a and ^ as seen at rof,. We 
place the observer at an arbitrary angle above the 
plane of Figure [2] Figure [5] shows the geometry 
and variables we will use to discuss the problem. 
The disk is represented by the circle. We define 
three coordinate systems. One corresponds to the 
disk (a;',y',z'), the z'-axis is normal to plane of 
the disk. The second reference (x, y, z) is centered 
at the disk but aligned in such a way that the 
observer is located on the x-axis at a distance Voh 
and y = y' . The point of observation Vob has an 
inclination l £ (0, 7r/2) with respect to the normal 
to the disk (see Figure [6]) . 




Fig. 5. — A neutrino is emitted from the point 
rem, which is at some distance above the plane of 
the disk (circle). The black hole is at the center 
of the disk. The x',y',z' correspond to the disk 
system, z' is normal to the disk. With respect 
to this system the emission point has coordinates 
(rem, 'fijf^ = T!'/'2 ~ (3)- V' IS thc angle between the 
x'-axis and a line joining the emission point to the 
black hole. 

The relation between these the coordinate 
frames is given by a rotation around the y-axis: 

Cx' = (sin t, 0, — cos i), (8) 
By' — e.y, (9) 
6*2' — (cost, 0, sin t). (10) 

In Figure[6]we show the third coordinate system 
which corresponds to the observer reference frame 
and is just a translation of the x,y,z reference with 
origin at rob- In this case Cdir = -x, Bright = y, 
&up = z. Note that Cdir points towards the black 
hole. 

An emitted neutrino has coordinates (rem, </5, ^ = 
7r/2 — /3) with respect to the x', j/', z' axes (Figure 
[5]). The angle between the x'-axis and rem is ip- 
On the other hand, on the x, y, z axes, the x-axis 
and the emission point form an angle ipem (Fig- 
ure |6|. The neutrino trajectories are restricted to 
a plane, called the observational plane, because 
of the spherical symmetry of space time in the 
Schwarzschild metric. This plane is defined by the 
null geodesic traveling from the emission point to 
the observer and the x-axis. We therefore can 
solve equation |3] on that plane, which corresponds 
to solving for a and ^ with initial conditions rem 
and (t>em- 

An image plane (as seen from the observer's ref- 
erence, i.e. consisting of the xy plane in Figures 
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Fig. 6. — Same coordinate systems as in Figure 
[5j The angle between the axis normal to the disk, 
z', and the observer direction x, is l. Here the 
coordinates of the emission point are given as seen 
from the x, y, z system. The neutrino trajectory 
forms a plane with the x-axis. The neutrino hits 
the observer plane and forms and angle a with 
respect to the Bright direction and an angle ^ with 
respect to edir- 

[5] and [6]) will intersect the observational plane at 
an angle a with respect to Crighti or equivalently 
the y-axis. The angle a can be obtained from the 
relationships between spherical triangles of Figure 
[7] (top). In this fi gure, the neutrino is emitted at 
fern and follows a trajectory represented by the 
dotted line reaching the observer at Vob- The an- 
gle between the neutrino observational plane and 
the xz plane is f — a. Using the spherical trian- 
gle defined by the arcs tp, 4> and f — i, and the 
spherical triangle defined by the arcs '0, 4> and (3 
we find 

sin(^cos/3 

Cosa= : — . (11) 

sm0 

Transforming the coordinate x' of the emission 
point to the x axis, after a rotation around the 
y axis by the angle t, we find 

COS0 = sin i cos -t- COS i sin /3 (12) 
sm(j) = \/l — cos^ (13) 

where cos '0 = cos /3 cos ip. 

Solving for the different cases of equation |3] we 
find the unique geodesic that connects the point 
{rem,V>em = to {rob,(t>ob = 0). This means we 
find h and a and by virtue of equation|4j we find ^. 
We therefore have a relation between (rem, f) and 
(^, a). Once ^ and a are found for every point over 



z' z 




Fig. 7. — Same coordinate systems as in Figure [5] 
Here we have changed the point of view in such a 
way that the observer's system seems rotated with 
respect to the disk system. The blue dashed circle 
represents the equatorial plane of the disk. The 
angles (/?, (/), -0, a and t are as defined in Figures [5] 
and[6) A neutrino is emitted at the rem and arrives 
to the point Voh- The neutrino trajectory (green 
dotted line) belongs to the observational plane. 
The observational plane forms an angle a with the 
y axis in the yz plane. Top: a relationship for is 
found by using the spherical triangles described by 
the arcs f/', and f — i; and /3, ^ and ^/i. Bottom: 
the neutrino trajectory forms and angle r\ with the 
plane x'y' . Using the spherical triangle described 
by the observational plane, the x' y' and the x' z' 
planes we find a relationship for r] in terms of a 
and i. 

the neutrino surface, we can calculate the solid 
angle dilob = sin^dfda covered by the observer. 

It is also useful to find a relationship for the 
angle 77 that the neutrino trajectory makes with 
the plane of the disk in terms of the already known 
angles. We find r] by using the law of cosines for 
the spherical triangle described by the neutrino 
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trajectory, the disk and the xz planes (see Figure 
[7] (bottom)) 



cos r] 



cos a sin t. 



(14) 



2.2. Energy red-shift 



The energy measured by an observer Eob 
changes from the emitted energy Eem by the red- 
shift factor (1 + z), 



E. 



1 + z. 



ob 



(15) 



In the case of a massless neutrino its mea- 
sured energy is the projection of the neutrino 4- 
momentum p on the 4-velocity u of the emitting 
matter, ppu^ . Then a general expression for the 
red-shift is, 

l + z=- T -2^. (16) 

(ptU' + PrW + P0U" + P^WP) 



In equation 16 (and in the next discussions) the 



subindexes em and oh indicate that the compo- 
nents of the metric tensor (or any other quantity) 
should be computed with the observer or emitter 
coordinates accordingly. From this general expres- 
sion we can consider specific cases according to the 
movement of the emitting particle and the metric 
of the space-time. For example, if the emitter is 
considered at rest then u*^ = = w'' = 0. For a 
cloud of nonrotating gas accreted by a black hole 
we have u'^ — while ^ 0. We could also 

have a case where the emitting and the receiving 
particles do not have movement in the r and d di- 
rections, and therefore — — '^i but 7^ 0. 
In what follows we study different cases depending 
on the emitter motion and the space-time metric. 

2.2.1. Non-rotating black hole 

For a non-rotating black hole we consider here 
two cases: in the first one the emitter and the 
observer do not have relative motion, and in the 
second case both rotate around the black hole. In 
the first case m'^ = = w'' = 0. Then the general 



expression of equation 16 becomes 



l + z 



pl"h 



pob^t 



(17) 



ob 



where pt is the neutrino energy which is constant 
along the neutrino trajectories. Therefore we have 



1 



em 

Kb ' 



(18) 



where m* = dt/dr. dr is calculated according to 
the metric describing the curvature of the space- 
time, 

,1/2 



dr 



-gajsdx^dx^^ 



(19) 



If the space-time around the black hole can be 
described by the Schwarzschild metric (equation 
[1]) then there is spherical symmetry, and therefore 
dr = {—gttY/'^dt. We then find the well known 
expression for the energy red-shift. 



l + z^ 



,1/2 
9tt)ob 



(-fftt) 



\l/2' 
[~gtt)em 



(20) 



where gtt — —(1 — rs/r). 

In the second case we are considering here, the 
emitter and the observer are rotating in stationary 
orbits around the black hole with angular veloci- 
ties ftem and ftob, respectively. This case exempli- 
fies a particle that has a stationary orbit around 
the black hole and is located at a fixed angle Oem- 
This particular choice could also represent ele- 
ments of matter of a steady state accretion disk. 
In Keplerian rotation the values for the angular 
velocity for the observer or emitter fl is 

Afi/2 

(21) 

In this case we have for the radial and azimuthal 
velocities = = 0, and therefore the gen- 



eral expression for the red-shift (equation 16 1 is 
reduced to 



l + z 



(22) 



The second fraction in the parenthesis equals the 
angular velocity = u'^/u* = dip/dt. Further- 
more, p^p and Pt are constants along the neutrino 
trajectories, which correspond to the projection of 
the angular momentum on the z'-axis Lz' and en- 
ergy E respectively. In the previous section we 
introduced the impact parameter b = L/E. Using 
these qviantities we re-write the red-shift as 



l + z 



(1 + ilenibcOSr]) 
M*^ (1 -I- Uobb COST]) 



(23) 
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where rj is the angle formed by and the plane 
of the disk 6* = 7r/2 (see Figure (I]). 




Fig. 8. — A neutrino is emitted with momentum 
P. The neutrino trajectory and momentum form 
an angle i] with the x'y' plane. The z' component 
of the angular momentum is L^' ~ Lcosi]. 

Using again the fact that u* — dt/dr, we get an 
expression for the energy shift in the Schwarzschild 
metric, when both observer and emitter rotate 
around the black hole (this derivation can be found 
in ( |Luminet||l979| or ( |Thorne|p7T| ), 



1 + z 



Clernb cos 7]) 



1 /2 

{~gtt)em (1 + ^obbcOS-q) 



(24) 



2.2.2. Rotating black hole 



In the case of a non-charged, rotating black hole 
the curvature of the line element can be written in 
the Kerr geometry as, ( Misner et al.|[l973 ), 

ds^ = gttdt^ +2gt^dtdip+g^^dip^ +grrdr'^ +gggde'^ , 

(25) 

with 



gtt 

9tv 
gee 



= - I- 



Y2 



tVa, 



r.ra sm 



sm 



f2 
^2 



(26a) 
(26b) 
(26c) 

- Aa^ sin^ 6'](^6d) 
(26e) 



Here a = Jc/GM'^ is the spin of the black hole 
(J is the total angular momentum), and A and T 
are given by 



A : 

'Y^2 



Tgr + a 
- a? cos^ I 



(27a) 
(27b) 



Considering the first case in which the emitter 
and observer have no relative motion we have a 
similar expression for the the red-shift as in equa- 
tion [20] but in this case with gu defined as in equa- 
tion [26al 



,1/2 

gu)ob 



{-gu) 



1/2 



>ob 



,1/2 
gttjem 



{-gtt) 



l- 



,1/2- 



(28) 



Now, as in the previous subsection, we consider 
the case where the observer and the emitter rotate 
around the black hole. The conditions for the 4- 
velocities {u^ = = 0) and the metric of space- 
time lead to 



dT 



dt^ 



-gtt - 2g, 



t^p 



dip 
dt 



gtpv 



d£_ 
dt^ 



1/2 



or in terms of the angular velocity, 

dt 1 
d7 



{-gtt - '2gt<p^ - g^^"^) 

Therefore the red-shift becomes, 



2^l/2■ 



(29) 



(30) 



1 



[-gtt - "^gtv^ 



ob 



02 1 1/^ 
yvV^^obiob 



[—gtt — "^gtip^em 



gipip^ernlJm 



(1 + Vtemhcosrf) 
(1 + Vtobhcosif) ' 



(31) 



Finally, if we assume that the emitter and ob- 
server are in Keplerian rotation, then their angular 
velocities are given by 



n = 



Ml/2 



'3/2 + aM^/'^ '' 



(32) 



with r evaluated consistently with the observer 
and emitter coordinates. The above analysis re- 
duces to the results of the Schwarzschild metric 
when we take the spin of the black hole a = 0. 

In Figure|9]we compare the energy shifts result- 
ing from taking different values for the spin of the 
black hole with a mass of 2.5M0 and allowing (or 
not) the emitter and the observer to orbit around 
it. In this example, we consider both observer and 
emitter located in the same plane xz {^p = Q). We 
also examine the effect of rotation by comparing 
two cases. We place emitters at the same distance 
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rem but located at two different different angles 
if = 90°, 270°. In can be seen from this figure 
that there is not a significant difference between 
using a Kerr or a Schwarzschild metric when cal- 
culating the energy shifts if the spin of the black 
hole is moderate. A larger effect, however, comes 
from taking into account the possible rotation of 
the observer and emitter around the black hole. 




Fig. 9. — Energy shifts for different black hole 
spin and rotation values compared to the non- 
relativistic case (No GR). The observer is located 
at (rob = 64.03 km, ipot = 0, 6* = 52.4°). The emit- 
ter is at (rem = 22.75 km, 6 = 48.75°). Unless 
indicated otherwise, ip^m — 0. Q — implies that 
the observer and emitter do not rotate around the 
black hole. The mass of the black hole is 2.5M0. 
The significant effect in the energy shift is due to 
the rotation of the emitter and observer around 
the black hole. 



2.3. Fluxes 

We are interested in determining the neutrino 
fluxes observed at some point in the vicinity of 
a black hole. We start by determining the emit- 
ted fluxes from the neutrino surface. We assume 
a Fermi- Dirac distribution which is suitable for 
neutrinos. At the emission point this is 



(33) 



with (7i^ = 1, Te„j the temperature at the emission 
point, and with an assumed neutrino chemical po- 
tential = 0. Using the known expression for 



the energy shifts, we can write this flux in terms 
of the energy measured by the observer 



(pemiEem = Eob{l + z)) = 

9uc El,{l + zY 

2TT^{hcf eMEobil + z)lTera) + 1 



(34) 



The temperature measured by the observer re- 
lates to the temperature at the emission point by 
Tem = (1 + z)Tob- If we knew Tob we could replace 
this expression and then the shifts factors would 
cancel out in the exponential of equation |34] to get 



where 

4>ob{Eob 



(t>ob{Eob) 



9vC 



.{Ee 



El 



1' 



(35) 



(36) 



2^^{hcfeMEob/Tob) 

is the observed flux. 

In general, the observed effective flux coming 
from a finite source is 



1 

47r 



dVLob X (t)ob{Eob), 



(37) 



where dftob is the solid angle that the source sub- 
tends as seen by the observer. In this case one can 
use (j)obs from Eq. [35| 

Another way to see this is to start from the 



quantity I/E^ which is an invariant (Misner et al 



1973j), where / is the specific intensity, therefore 



we can calculate the observed flux in terms of the 
emitted flux as 



1 

in 
1 

47r 



dflob X 



dflob X 



Iob{Eob) 



E. 



ob 



^em (-^em) 

{l + zfEob 



where I^m has the form. 



em 



exp{Eem/Te.. 



V 



(38) 



(39) 



because we have assumed a Fernii-Dirac distribu- 
tion. 

We would like to flnd reaction rates of neutrinos 
with matter around a black hole accretion disk. 
For this reason we want to calculate fluxes with 
energies in the observed system but using the val- 
ues of emitted temperatures that come from the 
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disk models. Therefore we make a change of vari- 
ables back to the observed energy to finally get 



1 

47r 



El, 



exp(£;„;,(i + z)/re„o + i' 



(40) 

for the observed effective flux. 

In order to determine dTlob we assign a coordi- 
nate frame centered at the observer, located at a 
fixed point Tob- Then dilob = sin^d^da, where ^ 
and a the angles discussed in section |2.1| and can 
be obtained by constructing null geodesies. Note 
that these variables (a, ^) are different from the 
ones describing a coordinate system centered at 
the black hole. 

In Figure [lO] we compare the fiuxes obtained 
when different general relativistic effects are con- 
sidered. The plot shows electron neutrino fluxes 
from a thin flat disk (adapted from a 3D hydro- 



dynamical model, as was used in Surman et al 



(2008)) and observed at a point with spherical co- 
ordinates (rob = 64.03 km, ^pob = 0, 6* = 52.4°). 
When only energy shifts are considered and we ig- 
nore rotation around the black hole, we find that 
the resulting neutrino fluxes are lower than the 
fluxes obtained for a non-relativistic case. How- 
ever, taking into account the bending of neu- 
trino trajectories makes the neutrino fluxes larger 
compared to the non-relativistic case for energies 
around 10 MeV, and larger compared to the case 
when the ray bending is ignored but the energy 
shifts are included. On the other hand, the high 
energy talk of the distribution is reduced. Adding 
the effect of rotation to the energy shifts decreases 
the flux at this spatial location in the model. How- 
ever, the effect of rotation can go in either direc- 
tion, i.e. at some points it increases the flux and 
at some points it decreases the flux. 

3. Calculations of Nucleosynthesis from 
Black Hole Accretion Disks 

We apply the general relativistic corrections to 
emitted neutrinos, as described in the previous 
section, from two specific models of black hole ac- 
cretion disks. The general strategy is to take cal- 
culated neutrino surfaces from the two different 
accretion disk models, and use these to find fluxes 
for points above the disk. Using these fluxes, we 
then find the nuclear products resulting from the 
interaction of the neutrino fiuxes with material 




E [MeV] 



Fig. 10. — Observed electron neutrino fluxes com- 
ing from a fiat disk and observed at a point lo- 
cated in spherical coordinates {rob — 64.03 km, 
ipob = 0, = 52.4°). The mass of the black hole 
is 2.5M0. The solid line corresponds to the case 
where general relativistic effects are ignored and 
where the angular velocity of the disk 57 = 
(No GR). The thin dashed line (NB) represents 
the fiux when only the energy shifts were consid- 
ered (fi = 0). The bold dashed line (B) includes 
the energy shifts and the bending of neutrino tra- 
jectories. The dotted line includes all the effects 
plus the rotation of the observer and the emitter 
around the black hole. 



out-fiowing the disk. 

We consider one dynamical and one steady 
state disk model. The first is based on a 3D hy- 
drodynamical model of a black hole and a neutron 



star merger studied by Ruffert and Janka (Seti- 
awan et al.|[2004{ [Janka et aLl[T999l [Ruffert fc 



Janka 2001). The second is based on a one di- 



mensional steady state disk modeled by [Chen &: 



Beloborodov (2007). 



The simulation by Ruffert and Janka corre- 
sponds to the merger of a I.6M0 neutron star and 
a 2.5M0 black hole with spin parameter a = 0.6. 
In this model general relativistic effects are in- 
cluded by using a modified Newtonian potential. 
The black hole is treated as a gravitational center 
surrounded by a vacuum sphere. The gravitational 
potential ^> bh of the black hole is an extension of 



the Paczynski-Wiita potential (Paczynsky & Wi- 



ita 1980 ) to a rotating black hole ( Artemova 1996 ). 
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As function of radius r, ^bh has the form 



BH 



GM 



BH 



dr 



(41) 



where /3 depends on the black hole spin parameter 
a, and ru, Mbh are the event horizon and mass 
of the black hole respectively. As a result of the 
coalescence, a disk is formed with inner boundary 
located at p = 14 km and surface extending to p = 
300 km (in cylindrical coordinates). The model is 
dynamical and therefore the disk we consider is 
based on a snapshot of the simulation. 



The model of Chen & Beloborodov (2007) cor- 



responds to a steady state disk. The mass of the 
black hole is SMq, the accretion rate M = 5Mq/s, 
and the spin parameters used are a = and 
a = 0.95. This model is fully relativistic. The 
disk is one dimensional, axially symmetric and is 
described by vertically averaged quantities. The 
disk extension goes as far as p = 600 km. For the 
vertical structure of the disk we used a simple hy- 
drostatic model that assumes an equilibrium with 
the gas radiation pressure and gravity. 

In Caballero et al. (2009) we calculated neu- 



trino surfaces in the corresponding 3D cylindrical 
grids for both the steady state and the hydrody- 
namical models. Our results for neutrino surfaces 
for a fixed angle (p can be seen in Figure 1 of |Ca-| 
ballero et al. (2009). The neutrino surfaces found 



are not smooth. They present sharp variations in 
temperature, density and height. This fact can 
be seen in the 3D image of the electron antineu- 
trino surface showed in Figure 3 of |Caballero et alT] 
(2009). We replicate here this figure for illustra- 
tion purposes (see Figure 11 ). On the other hand 
the neutrino surfaces from the steady state Chcn- 
Belobodorov model are smooth showing a torus 
shape. 

When determining the trajectories of the neu- 
trinos emitted from the disk, we consider two sce- 
narios. The first is the "flat disk" approximation. 
In this approximation, while the neutrino temper- 
atures are determined by the surface of last scat- 
tering, e.g. Figure [TT] the neutrino trajectories 
are started from the midplane (-z'=0) of the disk. 
The second is the "puffy disk" where the neutrino 
trajectories begin at the surface of last scattering. 

Since we are interested in the impact of the neu- 
trino general relativistic corrections on nucleosyn- 
thesis, we must consider the path of ejected ma- 



Fig. 11.— Figure 3 in (Caballero et al. 20091 



Electron antineutrino surface seen at some incli- 
nation angle (see the x, y, z axis on the lower 
left corner) . The height corresponds to the decou- 
pling height z^. The color scale corresponds to the 
emitted neutrino temperatures. The black area in 
the center represents the boundary with the black 
hole, r = 2rs. The neutrino surface shows a very 
uneven structure. 



terial from the disk. Here we consider two types 
of outflow trajectories. The first type is obtained 



with a similar prescription to that given in Sur- 



man & McLaughlin (20051, and used in Surman 



et al. (2008). In Figure 12 we show this outflow 



trajectory and the electron neutrino surfaces for 
the steady state and hydrodynamical models. The 
image corresponds to a transversal cut at an an- 
gle (fi = 0° . The trajectory starts at the point 
{p = 40, z = 48) km and extends vertically up to 
a turnover point at (p = 40, z = 96) km, after 
which it follows the radial direction. The outflow 
is taken to be adiabatic, with the velocity v of the 
outflow as a function of distance from the black 
hole r given hy v = Vad^ — fo/f)^, where tq is the 
starting position on the disk, v^o is the final coast- 
ing velocity of 0.1c, and fi determines the outflow 
acceleration, with lower /3 corresponding to more 
rapidly accelerating outflows. The second type of 
outflow trajectory considered is a parameterized 
spherical neutrino-driven wind trajectory, such as 



described in Panov & Janka (2009). This trajec- 
tory is purely radial, starting at {p = 40, z = 96) 
km. The outflow is again taken to be adiabatic. 
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but with a homologous velocity-radius dependence 
{v cx r) such that r(t) = roe*/"^, where r is the 
dynamical timescale. For both types of outflow 
trajectories, steady state conditions (i.e., r^pv = 
constant, where p is the baryon density) are as- 
sumed. More details of each outflow model can be 



(2009) 



found in Surman et al. ( 2008 1 and Panov & Janka 



N>300 



CB model v 




Fig. 12. — A transversal view of the neutrino 
surfaces at = 0° for the 3D hydrodynamical 
model(3D) and the steady state model (CB). The 
line OFl shows the trajectory followed by the out- 
flow. It starts in a straight line in the z direction 
and then at the turnover point follow the radial 
direction. The black ellipse a,t {x ~ y ~ 0) shows 
the black hole boundary. 

We calculate the observed neutrino fluxes at 
every point of the outflow trajectories, coming 
from all the points of the neutrino surfaces showed 



in Figure 11 and the corresponding torus-shaped 
neutrino surfaces of the steady state model. In or- 
der to calculate the fluxes we follow the neutrino 
trajectories emitted from the coordinates (p, ip, z^) 
(which can be translated to the spherical coordi- 
nates {Term fern, P — 7r/2~ 9 em) disCUSSCd in SCC- 

tion |2.1[ ), and arriving to the points of the outflow 
trajectory. Each point in the outflow trajectory 
has assigned coordinates r^b, fob, ^ob- In this way 
we can calculate the angles ^ and a, formed in 
the sky of every point of the outflow, by neutrinos 
traveling from the neutrino surfaces, and therefore 
we determine the solid angle described by the disk 
at those points. By means of Eq. |40] we calculate 
the corresponding fluxes. The energy shifts vary 



according to different conditions on the motion of 
the disk and the outflow, as well as the spin of the 
black hole. 

Using these neutrino fluxes we calculate the el- 
ement synthesis using a nuclear statistical equilib- 
rium code, a charged particle reaction code, and, if 
necessary, an r-process network code, as described 
in Surman et al. (20061. For the outflows from 
the asymmetrical 3D Ruffert and Janka disk, we 
take the additional step of following the outflow 
from starting points at four equally-spaced angles 
around the disk and then averaging the resulting 
abundance patterns. We consider a number of 
cases to assess the impact on the abundance pat- 
tern of the general relativistic corrections to the 
neutrino flux. 

4. Results 

As material flows away from the disk, it begins 
as dissociated free neutrons and protons. Then 
as the material cools, heavy nuclei are formed. If 
the conditions are right then very heavy elements, 
such as the r-process elements or p-process ele- 
ments, are formed. Of key importance in deter- 
mining the type of nuclei formed are the relative 
numbers of neutrons and protons. Neutrons are 
converted to protons through the weak interac- 
tion, in particular. 



Ve + n -ir^ p + e 



(42) 
(43) 



We consider matter outflows that begin at the 
surface of a trapped source of neutrinos and end 
far from the source. Close to the surface the neu- 
tron to proton ratio is determined by all four of 
these reactions. However as the material flows 
away, there is a period where the most important 
reactions are the electron neutrino and electron 
antineutrino capture reactions. 

In the top panel of Figure[T3j one can see this ef- 
fect in matter outflows from a black hole accretion 
disk. The solid lines show the electron fraction 
Ye = TT-pl (rip -|- n„), where Up and n„ are the num- 
ber densities of the neutrons and protons in the 
material. The dashed lines show the equilibrium 
electron fraction in the presence of only neutri- 
nos, i.e. what the electron fraction would be if (1) 
the neutrino and antineutrino capture reactions 
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Fig. 13. — The top panel shows electron frac- 
tions (solid lines) and neutrino equilibrium elec- 
tron fractions (dashed lines) for the spherical wind 
trajectory with an entropy per baryon s/k — 20 
and a timescale of r = 5 ms. In all cases the "flat 
disk" approximation to the black hole neutron star 
merger model was used so that all the neutrinos 
are launched from the z=0 plane. The yellow lines 
treat the neutrinos as Newtonian. The green lines 
correspond to a situation where general relativis- 
tic neutrino redshift and trajectory bending are 
included, but rotation and black hole spin are not 
taken into account. The red lines add the effect of 
disk rotation. The blue lines include a black hole 
spin of a = 0.6 into the Kerr metric when calcu- 
lating the gravitational redshift. The solid lines in 
the bottom panel show final nuclear abundances 
for each calculation, while the crosses show the 
scaled solar data. 

were strong enough to establish an equilibrium be- 
tween neutrons and protons and (2) electron and 
positron capture can be neglected. The different 
color lines correspond to different calculations of 
the neutrinos using the same disk model and mat- 
ter outflow. (The differences in the calculations 



are explained below.) It can be seen in the case 
of the yellow lines that the actual electron frac- 
tion approaches the equilibrium electron fraction 
at large distances. In the case of the red, green and 
blue lines, the neutrino and antineutrino fluxes are 
lower and this equilibrium is not attained. 

One can see also in this figure that the neu- 
trino equilibrium electron fraction is not always 
the same. This is due to the relative difference in 
the spectra of the neutrinos and the antineutrinos. 
If the antineutrinos have higher energy than the 
neutrinos, then the material will be driven neutron 
rich, and vice versa. 

In Figure [13] we start with the black hole neu- 
tron star merger model and examine the simplest 
case: a flat disk. When calculating the neutrino 
trajectories, we take the case of a black hole with- 
out spin, a = 0. The choice of a = leads us to 
calculate the neutrino trajectories from the disk to 
the points of the outflow using the Schwarzschild 
metric as we discussed in section 12.11 We find 
the energy redshifts using Eq 



20 of section 2.2 



The impact of the neutrino gravitational redshift 
and neutrino trajectory bending can be seen in 
the figure. For these calculations we have used a 
spherical neutrino-driven wind trajectory with a 
low entropy per baryon s/k = 20 and fast outflow 
(t = 5 ms). The yellow line shows the case with 
no neutrino general relativistic effects. The neu- 
trino and antineutrino fluxes are strong enough so 
that neutrino equilibrium of the electron fraction 
is obtained. Furthermore, the antineutrino flux 
has higher energy than the neutrino flux and so the 
material remains neutron rich and the first peak of 
the r-process abundance pattern is obtained (yel- 
low line bottom panel). 

In the same figure, the green line shows a calcu- 
lation where neutrino trajectory bending and neu- 
trino redshift have been taken into account but 
without including the effects of the rotation of the 
disk. One can see that the solid and dashed lines 
never meet, and the equilibrium neutrino electron 
fraction is never obtained. This is because both 
neutrino and antineutrino fluxes are less energetic 
when general relativistic effects are included. In 
particular the high energy tails of the neutrino and 
antineutrino spectral distributions are reduced. 

One can also see that the equilibrium electron 
fraction (green dashed line) is significantly higher 
than in the Newtonian case. This happens because 
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the antineutrinos are on average emitted closer to 
the black hole than the neutrinos. Thus the an- 
tineutrinos are more redshifted, and the balance 
of the weak rates moves toward proton rich ma- 
terial. However, in this fast outflow scenario, the 
first effect dominates. The equilibrium electron 
fraction is never obtained due to the overall de- 
creased strength of the neutrino and antineutrino 
fluxes. Therefore the material retains much of its 
original neutron richness, and the r-process pro- 
ceeds a little further - out to the second peak in 
this scenario (green line of the bottom panel). 

The remaining two lines in Figure [13] examine 
two additional effects. The red lines include the 
effect of disk rotation in the neutrino energy shifts, 
and the blue lines show the influence of the black 
hole spin (a = 0.6) on the redshift of the neu- 
trinos (the Kerr metric was not used to compute 
null geodesies). As can be seen from the figure, 
these effects create a more modest impact on the 
final abundance yields than the difference between 
treating the neutrinos as Newtonian or using the 
Schwarzschild metric to describe their evolution. 

In Figure [T3j much of the behavior of the abun- 
dance pattern shown in the bottom panel is due 
to the low entropy and fast outflow of the mat- 
ter. In Figure 14 we show the same model with 



the one exception that we replace the matter out- 
flow trajectory with one that has a higher entropy, 
s/k = 75, and slower outflow, t = 50 ms. It can 
be seen that the initial electron fraction is much 
higher due to the higher entropy. The final elec- 
tron fraction in the Newtonian case is similar to 
the previous example since the neutrino and an- 
tineutrino fluxes are high enough for equilibrium 
to be reached. In the cases where general relativis- 
tic effects are considered, again the actual electron 
fraction never obtains the equilibrium value. How- 
ever, in this case, the electron fraction remains 
higher than the equilibrium value. While these 
high entropy matter outflows produce interesting 
nucleosynthesis, an r-process is not produced. 

When considering neutrino general relativistic 
effects, the emission point of the neutrinos is cru- 
cial. Previous nucleosynthesis calculations from 
accretion disks have considered all neutrinos as if 
they were emitted from the z' = plane of the 
disk, i.e. a flat disk approximation. In Figure 15 
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Fig. 14. — The top panel shows electron frac- 
tions (solid lines) and neutrino equilibrium elec- 
tron fractions (dashed lines) for the spherical wind 
trajectory with an entropy per baryon s/k = 75 
and a timescale of r = 50 ms. In all cases the 
"flat disk" approximation was used so that all the 
neutrinos are launched from the z=0 plane. The 
line colors represent the same cases as in Figure 

m 



we show the consequences of abandoning the "flat 
disk" approximation. The same model and matter 



outflow trajectory was used as in Figure [13] but 
the "flat disk" approximation was not used and 
neutrino trajectories were started from the neu- 
trino decoupling surface. 

As can be seen from Figure |15[ the behav- 
ior trend seen in the "flat disk" approximation 
has been accentuated in the "puffy disk" scenario. 
The neutrinos are emitted further from the black 
hole and thus redshift effects are smaller, thus 
one might have expected the opposite behavior. 
However, the hottest neutrinos are emitted on a 
surface that tilts toward the black hole, and so 
the difference in flux from geometrical consider- 
ations more than compensates for a softening of 
the redshift. We note also that the electron neu- 
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Fig. 15. — Same as Figure 13 except that the "flat 
disk" approximation was not used. 



Fig. 16. — Same as Figure 14 except that the "fiat 
disk" approximation was not used. 



trino surface is higher than the electron antincu- 
trino surface. This enhances the flux of neutrinos 
relative to antineutrinos, and contributes to the 
proton rich equilibrium electron fractions seen in 
the top panel of Figure [15] Nevertheless, due to 
the overall reduction in strength of the neutrino 
and antineutrino fluxes, the material remains even 
more neutron-rich and a more robust r-process is 
obtained. 

If we instead compare the flat disk approxima- 
tion, Figure[T4j to the full "puffy disk" calculation. 
Figure |16[ using high entropy and slower outflow 
conditions, we find that neutrino general relativis- 
tic effects can cause material that would be neu- 
tron rich in a Newtonian calculation to become 
proton rich. 

To better understand whether these results are 
particular to a given disk and outflow model, we 
compare calculations done with the black hole neu- 
tron star merger model to those done with the 
steady state disk model. We use the parameter- 
ized wind with the full outflow trajectory shown 



in Figure 12 a low entropy per baryon s/k = 10, 
and acceleration parameter /3 — 0.2. In Figure [T7| 
we show the abundance patterns obtained in this 
way. The qualitative pattern seen in the results 
is similar to that seen in the low entropy cases of 
Figures [13] and [T5j 

5. Conclusions 

We have studied the influence of neutrino gen- 
eral relativistic effects on the spectra of neutri- 
nos after they are emitted from black hole accre- 
tion disks. We examined two models, a snapshot 
of a three dimensional dynamical calculation of a 
black hole neutron star merger and a steady state 
one dimensional disk. We find that both redshift 
and trajectory bending are significant, and will 
influence not only the spectra of neutrinos but 
also nucleosynthesis products from material that is 
ejected from the disk near the source of neutrinos. 

The overall impact of general relativistic effects 
on the neutrino spectra as compared with a New- 
tonian calculation is to reduce the energy flux of 
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Fig. 17. — Compares the results for the two differ- 
ent disk models using the outflow trajectory shown 
in Figure 12 with an entropy per baryon s/k = 10 
and acceleration parameter /3 = 0.2. The colors 
indicate the same treatment of the neutrinos as in 
the preceding four figures, although for blue line 
in the bottom panel a black hole spin parameter of 
a=0.95 was used. The top panel shows the results 
of a calculation which used the black hole neutron 
star merger disk and the bottom panel from a cal- 
culation which used the steady state disk. The 
"puffy disk" was used in both panels. 

the neutrinos, and in particular to supress the high 
energy tails of the spectra. Since the electron an- 
tineutrinos are emitted closer to the black hole 
then the electron neutrinos, the effect is stronger 
in this channel. Thus, the electron antineutrino 
energy flux is reduced more than the electron neu- 
trino flux. 

Neutrinos are emitted from their surface of last 
scattering. The conditions at the point of last 
scattering initialize the neutrino spectral proper- 
ties during free streaming phase. From the point of 
view of determining the flux at every point above 
the disk (i.e. integrating over the emission sur- 
face) , it can be computationally convenient to col- 
lapse the surface to the z = plane of the disk. We 



compared two approaches, where neutrinos with 
the same spectral properties were emitted from 
either the plane of the disk ("flat disk" approx- 
imation), or from the actual decoupling surface 
("puffy disk"). We found that the difference on 
the neutrino spectra above the disk was consider- 
able, with the "puffy disk" scenario resulting in a 
lower energy flux that the "flat disk" scenario as 
calculated at the same point above the disk. In 
the "puffy disk" scenario, there is less redshifting, 
since the neutrinos are emitted farther from the 
black hole. However, the effect of the geometry of 
the disk is strong. The hottest neutrinos are emit- 
ted on a surface that is sloped toward the disk and 
the neutrinos bend toward the black hole, which 
reduces the fluxes of both neutrinos in directions 
away from the black hole, but has a larger impact 
on the antineutrinos. Further, the neutrino emis- 
sion surface is higher, so that the neutrinos are 
emitted closer to points of nucleosynthetic inter- 
est above the disk than the antineutrinos are. This 
reduces the relative flux of antineutrinos relative 
to neutrinos at many points above the disk. 

We considered the impact of these effects on the 
element synthesis in matter ejected from the disk 
in the vicinity of the neutrino decoupling surface. 
The largest impact is in the reduction of the neu- 
trino fluxes. The neutrinos lose some of the influ- 
ence that they had in setting the neutron to proton 
ratio. Thus, the amount of heating that material 
receives as it leaves the surface of the accretion 
disk becomes an important factor in determining 
the final abundance pattern. In our models, with 
even modest heating and fast outflow, the neutron 
to proton ratio will be set almost entirely by elec- 
trons and positions, and only in the case of little 
to no heating will an r-process be produced. In 
the case of slower outflow, the neutron to proton 
ratio will be set by the neutrino fluxes, but since 
the antineutrino flux is reduced more by general 
relativistic effects than the neutrino flux, the elec- 
tron fraction of the material ranges from 0.4 to 
0.6, which will produce interesting nucleosynthe- 
sis, but not an r-process. 

Our study suggests that a range of nucleosyn- 
thesis products may be possible for accretion disks 
with trapped electron neutrinos and antineutri- 
nos, and that the outcome is dependent on the 
neutrino spectra. Neutrino spectra are sensitive 
not only to the general relativistic effects outlined 
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here, but also to the decouphng surface of the neu- 
trinos. Therefore the nucleosynthetic outcome in 
ejecta from black hole accretion disks is as well. 
Thus, more detailed neutrino diffusion, preferably 
including the effects of general relativity, is war- 
ranted in future studies. 
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